NASA Technical Paper 1261 


; NASA 
' TP 

I 1261 
0,1 


copy. 

4^1, 

Procedures for Generation ’ ° ^ 

and Reduction of Linear Models 
of a Turbofan Engine 


■ I 


1=3 I 
b-* ! 
UJ i 


i UJl 

! x: j 

ii 03! 


Kurt Seldner and David S. Cwynar 


AUGUST 1978 




TECH LIBRARY KAPB, NM 



NASA Technical Paper 1261 


tech library kafb, nm 



□13M3MB 


Procedures for Generation 
and Reduction of Linear Models 
of a Turbofan Engine 


Kurt Seldner and David S. Cwynar 
Lewis Research Center 
Cleveland, Ohio 


fVASA 

National Aeronautics 
and Space Administration 

Scientific and Technical 
Information Office 


1978 


PROCEDURES FOR GENERATION AND REDUCTION OF LINEAR 
MODELS OF A TURBOFAN ENGINE 
by Kurt Seldner and David S. Cwynar 
Lewis Research Center 

SUMMARY 

A real-time hybrid simulation of the Pratt & Whitney FIOO-PW-FIOO turbofan en- 
gine was used for linear-model generation. The linear models were used to analyze 
the effect of disturbances about an operating point on the dynamic performance of the 
engine. A procedure that disturbs, samples, and records the state and control vari- 
ables was developed. 

For large systems, such as the FlOO engine, the state vector is large and may con- 
tain high-frequency information not required for control. Thus, reducing the full-state 
to a reduced-order model may be a practicable approach to simplifying the control de- 
sign. A reduction technique was developed to generate reduced-order models. Selected 
linear and nonlinear output responses to exhaust-nozzle area and main-burner -fuel -flow 
disturbances are presented for comparison. 

INTRODUCTION 


Durii^ the past decade linear- control-design techniques for multiple-input - 
lultiple- output (multivariable) systems have been the subject of much research. Most 
[■equire a linear time domain or state-space formulation of the dynamic characteristics 
3f the process. Although modern turbofan and turbojet engines are considered multi- 
^•ariable processes, they are normally described by many nonlinear, d3mamic relations. 
Thus, linearization procedures must be used to obtain a linear, state-space representa- 
jtion of the nonlinear, multivariable process. Because the linear models are only valid 
^bout a particular operating point, several models must be generated to describe the 
^ngine operation over a wide operating envelope. Once linear models are obtained, the 
control can be designed using a linear multivariable procedure. 

Digital and hybrid engine simulations are available to extract linear models at 
Various operatii^ points throughout the flight envelope. This report describes the pro- 



cedure for extracting linear models from a real-time hybrid computer simulation of a 
turbofan engine - the Pratt & Whitney FIOO-PW- 100(3) engine (ref. 1). The computer 
simulation was developed as part of the FlOO multivariable -control, synthesis program, 
sponsored jointly by the Air Force Aeropropulsion Laboratory and the NASA Lewis Re- 
search Center. 

The report includes (1) a description of the engine simulation, (2) a brief descrip- 
tion of the analytical basis for the linear-model generation technique, (3) a method for 
reducing the resultant linear-model complexity by eliminating unwanted information, 
and (4) a presentation of selected d 3 mamic responses for both nonlinear and linear 
models. The hybrid-computer listings for generating the disturbance data for comput- 
ing linear models, system matrices, and eigenvalues for several operating conditions 
are presented in the appendix. 


DESCRIPTION OF FlOO TURBOFAN ENGINE 

The Pratt & Whitney FIOO-PW- 100(3) engine (fig. 1) is an axial, mixed-flow, twin- 
spool, low-bypass-ratio turbofan with afterburning capability. The engine has a single 
inlet for both fan and core airflow. The fan airflow is separated into two streams: one 
passing through the fan duct and the other passii^ through the engine core. A three- 
stage fan is connected to a two-stage, low-pressure turbine. A 10-stage compressor 
is connected through a hollow shaft to a two-stage high-pressure turbine. The fan has 
variable guide vanes; the compressor has a variable guide vane followed by two vari- 
able rear stator vanes. Compressor discharge bleed air is used for cooling the high- 
and low-pressure turbine blades and for powering the augmentor turbopump. 

The main combustor consists of an aimular diffuser with 16 fuel nozzles. The duct 
and core streams combine in the augmentor section and discharge through a variable 
convergent-divergent nozzle. The augmentor consists of a diffuser section with five 
concentric fuel manifolds. 

The basic engine-control system consists of a hydromechanical fuel-control system 
and an electronic-supervisory-control system. The function of the hydromechanical 
control is to (1) meter fuel to the main burner as a function of power-lever angle (PLA), 
compressor speed Ng, fan-discharge total pressure P^g, and compressor-discharge- 
static pressure ^S3» (2) position the compressor-stator vanes to improve starting and 
high Mach number characteristics, (3) meter fuel to the five afterburning zones as a 
function of power -lever angle, fan-duct discharge temperature, and compressor- 
discharge- static pressure, and (4) control exhaust-nozzle area to maintain desired en- 
gine airflow during augmented operation. 

The electronic supervisory-control positions the inlet -guide vanes, trims main- 
combustor-fuel flow to maintain engine operation within safe limits, and trims exhaust- 
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nozzle area to satisfy engine-flow requirements (refs, 1 and 2). 


ENGINE SIMULATION 

The development and evaluation of controls for aircraft propulsion systems requires 
a dynamic model of the engine and control system. The simulation allows us to study 
and predict engine performance over the specified operating range. The mathematical 
model can also be used to evaluate various control concepts and to study problems dur- 
ing the development and testii^ phases of an engine. 

The details of a hybrid- computer simulation for the P&W FlOO- 100(3) turbofan en- 
gine are reported in reference 1. The engine was modeled in real time on the EAI 
680-640 hybrid computer. Real-time simulation permits operation with d^ital control 
software. Figure 2 is a block diagram of the interconnection of the various ei^ine com- 
ponents as well as their significant variables. The real-time, hybrid- computer, ermine 
simulation was modeled according to the engine manufacturer's digital simulation. The 
steady-state characteristics of the rotating elements (compressor, fan and turbine per- 
formance maps) were extracted from this digital simulation. The digital portion of the 
hybrid computer was used to generate nonlinear bivariate functions for the component 
performance maps. Interconnecting volumes between components were used whenever 
the gas dynamics are significant or are required to prevent computer instability due to 
algebraic loops. The continuity, energy, and state equations were solved for each 
volume. The equations for conservation of angular momentiun are used for the fan- 
and compressor-speed calculations. 

LINEAR MODEL GENERATION 
Linearization of Model 

The FlOO engine is a complex, nonlinear, dynamic process consisting of many 
states. Small-signal-linearization techniques (ref. 3) can be used to represent the non- 
linear process about an operating point by a set of first-order matrix differential equa- 

( 1 ) 

( 2 ) 

(2a) 
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tions. The system equations can be written as 


X = Ax + Bu 


y = Cx + Du 
z = Hx 


Illlllllllllllll 


Equation (1) defines the relation between the system states (i = 1, 2, . . . , n) and 
the forcing function Uj (j = 1, 2, . . • , p). Equation (2) relates the outputs y to the 
states Xj and control inputs Uj. The symbols are defined in appendix A. 

Linear models generated from the real-time, hybrid- computer, engine simulation 
contain all states defined by the engine model. The states as well as control and output 
variables are as follows: 


Engine state variables: 

(1) Fan speed, N^, rpm 

(2) Compressor speed, Njj, rpm 

(3) High-turbine “inlet pressure, N/cm^ 

(4) Stored mass in the fan- turbine -inlet volume, kg 

(5) Nozzle-inlet pressure, P7» N/cm^ 

(6) Fan -turbine -inlet pressure, P4P N/cm^ 

(7) Stored mass in the fan duct, Wjg, kg 

(8) Fan-duct temperature, T^g, K 

(9) Stored mass in the high-turbine-inlet voliune, W^, kg 

(10) Stored mass in the compressor-discharge volume, Wg, kg 

(11) Compressor -discharge temperature, T^, K 

(12) Afterburner-inlet pressure, Pg, N/cm^ 

(13) Stored mass in the afterburner -inlet volume, Wg, kg 

(14) Stored mass in the nozzle-inlet volume, W,^, kg 

(15) Afterburner gas flow, Wg, kg/sec 

(16) Fan-duct airflow, Wjg, kg/sec 


Control variables: 

(1) Fan-inlet-guide -vane angle, CIVV, deg 

(2) Rear compressor variable vanes angle, RCVV, deg 
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(3) Nozzle throat area, A^^, m 

(4) Main -burner -fuel flow, w^^, kg/hr 

(5) Afterburner fuel flow, w^.^, kg/hr 


Output variables: 

(1) Fan airflow, w^ ^^ , kg/sec 

(2) Rear-compressor-variable-vane angle, RCVV, deg 

(3) Net thrust, F^, kN 

(4) High-turbine-inlet temperature, T^, K 

(5) Fan-discharge-flow parameter, AP/Pj^g 

(6) Fan sail margin, SMF 

(7) High-compressor-stall margin, SMHC 
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Linear models must be generated about an operating point aloi^ the normal operating 
line of the ei^ine. Since these models only describe the engine performance for small 
input disturbances at a specified operating point, several linear models must be gener- 
ated. The process of disturbing, sampling, recording, and reducing the raw data is 
automated. Disturbances of state and control variables are used to generate data from 
which the various elements of the system matrices of equations (1) and (2) can be deter- 
mined. The magnitude of the disturbance is limited to ±2 percent to insure that per- 
formance is maintained within a linear range of operation. The resultii^ linear models 
describe the dynamic performance of the engine about the specified operating point. 


Linear Model Generation 

The process of acquirii^ the data for model generation is outlined in figure 3. To 
generate linear- state models, the engine simulation must be set up at a specified power 
and flight condition. At the desired operating condition a hybrid subroutine (STATE 1) 
can be used to set the initial conditions for the desired operating condition. After the 
initial conditions (steady- state values for state variables) are set, the anal(^ portion of 
the hybrid computer can now be placed in the ’’initial condition” mode. This procedure 
prevents integrator drift during the sampling and data recording process. Another hy- 
brid subroutine (STATE 2) applies the disturbances to the state and control variables. 

The disturbance data are trunked to the SEL810B digital computer for sampling and re- 
cording. The two subprogram listings (STATE 1 and STATE 2) are given in appendix B. 

The block diagram in the lower portion of figure 3 illustrates the sampling, verify- 
ing, and recording processes of the disturbance data. The processes are outlined in 
more detail in the flow charts of figure 4. On initiation of the sampling process, a time- 
delay circuit is activated to insure that the simulation is steady state. A latch circuit is 
also set to L = -2 and updated during the process for logical decisions. Each engine 
variable was sampled 4096 times, and a mean was calculated for each variable. 

The sampling process was repeated, and a second mean M 2 was calculated. The two 
means Mj^ and M 2 were compared. If the error between the means is within 5 milli- 
volts, the data are accepted and descaled to engineering units. If the error between M^ 
and M 2 exceeds 5 millivolts, the sampling and validation process is repeated. At this 
time the latch circuitry is also updated. If the Mj^ and M 2 are not within tolerance 
after the second trial, the process is terminated and the operator alerted. All accept- 
able data are stored on a floppy disk before being transmitted to the IBM 360 computer. 

The sampled data must be read into the IBM 360 computer. Figure 5 presents a 
block diagram showing the data verification, processing, and linear model generation. 
After the sampled, descaled data are read into the computer, the mantissa and charac- 
teristic of each variable are added to form a check-smn. The check-sum value is 
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compared with that stored on the floppy disk by the SEL810B computer. If the two 
values agree, the process continues to the validation mode. If the two values do not 
agree, the process is terminated, and the data must be corrected. 

The data were then validated by verifying that the correct variables were disturbed 
and that the removal of the disturbance resets the state and control variables. A base 
run (without disturbances) was used for comparison. The tolerance between the two 
sets of data (base and reset) must be within a specified tolerance. If they are within 
tolerance, the data can be processed to compute the system matrices. The matrix com- 
putations used both positive and negative disturbances to obtain a better average and 
minimize nvunerical errors. The system matrices were stored in the IBM 360 computer 
for later use in calculating the dynamic characteristics of the linear models. 


LINEAR MODEL OPERATIONS 
State Vector Transformation 

Some of the states generated by the FlOO, real-time simulation, such as mass, 
may not be measurable and must be converted to either a pressure or temperature re- 
lation. Taylor’s series and the equation of state were used to modify the state vector 
so that it included only desirable measurements for the linear models. These generated 
models are valid only for small disturbances about a steady-state operating point. The 
mathematical analysis relating the linear transformation of the state vector is as fol- 
lows: Expand Taylor's series about the nominal operating conditions Pq’ ''^0’ 
where these values are the steady-state conditions at the engine-operating point. The 
Taylor-series expansion for temperature T about the nominal steady-state point Tq 
becomes: 


(T - Tq) = -^ 
^ 3P 


(P 


P=Pn 


^ 8W 


(W 


P=Pr 


Wo) + — 

^ 2! 


8^T 

8P^ 



W=Wq 

|w=Wq 

1 

1 

.2 

(P - 

Po)(w-Wo)+-^t 

(W-WJ^ 

3P aw 

P=Pq 

W=Wq 

\J ey 

P=Pq 

W=Wq 



(P - Pn)' 


P=Pq 

w=w, 


0 


+ higher order terms (HOT)I 


(3 
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The equation of state is 


P = — WT 


(4) 


where is the gas constant of air (in N • cm/kg - K) and V is the volume (in cm'*). 
We may rewrite equation (4) so that 


T_ V P 
■r^w 


(4a) 


Differentiating equation (4a) yields 


0T 


dT 


dP R^w aw 


V P 
W^ 


a^T 


ap^ 


= 0 


a^T 


= 2 


a^T 

ap aw 


V p 


RaW' 


aw^ ^A w^ 


And substituting into equation (3) yields 


(T - Tq) = (P - Pq) - (W - Wq) 

U 1-* T17 U o W 


R^W 


R^W‘ 


-X— (p - Pq)(W - Wq) + (W - Wq)^ + HOT 


rt u V o 

R^W^ R^W-^ 


If we neglect second- and higher-order terms, equation (5) becomes 


(T - To) = (P - Pq) - (W - Wq) 


R^W 


R^W‘ 


If 


K = AT = (T - Tq) ap = (P - Pq) aw = (W - Wq) 


R 


(5) 


( 6 ) 


then, equation (6) reduces to 


AT = — AP - — AW 


W 




( 7 ) 


Differentiating equation (3) yields 


AT = -^ AP + AW + AP ■ AP + AP AW 


dP 


aw 


ap^ 


ap aw 


+ AP AW + AW AW + HOT 


ap aw 


( 8 ) 


dW 


When all second-order and higher-order terms are neglected, equation (8) reduces to 


AT = -^ AP + -^ AW 


ap 


(9) 


aw 


Substituting for the partial derivatives yields 


AT = — AP - — AW 
W 


( 10 ) 


Similarly, for AP 


ap 


ap 


AP = AT + — AW 

aT aw 


( 11 ) 


or 


AP = ^ AT + ^ AW 


K 


K 


( 12 ) 


Using equations (10) and (12) allows the state vector to be transformed to include only 
pressure, temperature, and flow relations. The transformation of the state vector can 
be achieved through a transformation matrix E that relates pressure and temperature 
to mass in the volume. In matrix form 
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(13) 


Xj = Ex 

where X is the original state vector, Xj is the transformed state vector, and E is 
the transformation matrix, Substitutii^ equation (13) into equations (1), (2), and (2a) 
and simplifyii^ yield 


X^ = AX^ + Bu (14) 

y = CXj + Du (15) 


z = HXj 


(16) 


where 


A = EAE~^ n X n 

B = EB n X p 
C = CE"^ rxn 

and where A, B, and C are the original state, control, and output matrices. 

Equations (14) to (16) are linear models composed of measurable states. The 
transformation is a similarity transformation so that the ordinal set of eigenvalues 
will be retained. The transformation matrix E, of course, varies with engine operat 
ii^ condition; H is the identity matrix since the state vector Xj contains the trans- 
formed states. 


Model-Order Reduction 


Before reducing the order of the linear model, it may be necessary to permute the 
state vector so as to group slowly and rapidly decaying states separately. Permutation 
of the state vector may also be used to select particular states for the lower order 
Imodels. A similarity transformation can be used to perform this procedure: Let 


*2 = Vl 


( 17 ) 


where Xj^ is the transformed state vector, X 2 is the permuted state vector, and Ep 
is the permutation matrix. Substituting equation (17) for state vector X^ in equations 
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(14) to (16) yields the permuted linear model: 


where 


^2 = ^ 
y = C„X2 + Du 


( 18 ) 

(19) 

( 20 ) 




®n ^ ^ P 


C„ = CE"^ r X n 
n p 

Hn = HE“^ n X n 


and where A, B, and C are the transformed state, control, and output matrices. 
Equations (18) to (20) form the linear model for a reordered state vector. If permuta- 
tion of the state vector is not desired, the resulting matrix E„ is the identity matrix. 

-1 -1 P 

The output matrices HEp and CEp transform the permuted state vector to the orig- 
inal state output and output vectors. This procedure is required to maintain the original 
order of the outputs for data identification and plotting. 

The selection of an appropriate set of states for the reduced model is the problem 
imposed by linear model reduction. Model-order reduction can be based on the princi- 
ple that the state vector can be permuted into separate groups of slowly and rapidly de- 
caying states. (Usually a 10: 1 separation between eigenvalues can be considered as ac- 
ceptable. ) However, the eigenvalues in most physical systems, caimot be readily 
separated, and the identification of states with specific eigenvalues is a difficult, if not 
impossible, task. In addition, specific states must be included in the reduced-order- 
linear-model representation. Therefore, the reduced-order-linear-model state vector 
must be based on intuition and turbofan- engine controls experience. The states that 
should be considered are those required by the engine controller and protective devices. 

The state vector must be partitioned into two state vectors: one with the desired 
states for the reduced-order model, and the second with the remaining states that will 
be eliminated. The dynamic characteristics of the latter states are neglected by setting 
the derivatives to zero. These states are approximated by their steady-state values. 


The eliminated states are then reconstructed from information contained in the reduced- 
order model. The reduction technique is described in references 4 to 6. The mathe- 
matical formulation is as follows: 


X 


s 



( 21 ) 


where Xgj is state vector of retained system states and X ^2 is the state vector of 
eliminated system states. Applying the notation of equation (21) to equation (1) yields 



( 22 ) 


where Ajj, Aj 2 > ^2V ^22* ®1* ®2 matrices resultii^ from the partitioning 

of the n X 1 state vector X^. The solution of equation (22) for X^j and Xg 2 is 


^sl " -^ll^sl '^12^s2 ®l" 

^s2 ^ ■^21^sl ^ ^22^s2 ^2^^ J 


(23) 


By setting X ^2 = ® substituting into the Xgj 


relation 


^s2 " "'^22^‘^21^sl ^ 2 '^^ 


(24) 


the Xgj relation becomes 


^sl “ (^11 “ '^12^22‘^2l)^sl (®1 ■ '^ 12 ^ 22 ^ 2 )'^ 


(^1 


For the output equation (2) 




s2y 


(25) 


(26) 


where Cj^ and C 2 are matrices resultii^ from the partitioning of the output vector y. 
Substituting for ^s2 in equation (26) yields 
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II 


Similarly, for equation (2 a) 


= (- 


H, - H2A‘^A2i)Xgj - H„A;iB„u 


^2^22^2' 


(28) 


where and H 2 are matrices resulting from the partitioning of the state output vec- 
tor z. The dimensions of the state, control, and output vectors, and system matrices 
are given in appendix A. Rewriting equations (25), (27), and (28) yields 


" Vsl * V 




(29) 


where 


^ = ^s^sl V 


“^s = (^11 ■ ■^12'^22^2l) 

®s ^ (®1 ■ ■^12-^22^2) 

^s ^ (^1 ■ ^2‘^22‘^2l) 

^s = (d - C 2A22S) 

^s "" (**1 ■ ®2^22'^2l) 


Hu = -H 2 A 22 H 2 


Equations (29) represent the reduced-order-linear model. The system matrices 
(A, B, C, D, H) must be replaced by the reduced-order -linear model matrices (A^, Bg,| 
^s» ^s» Hg). The reduced-order model matrices contain sufficient information to re- 
construct all eliminated states and outputs. It should be emphasized that the reduced- 
model accuracy is strongly dependent on the eigenvalues for the retained and eliminated 
states. If the ratio of the lowest eigenvalue for the eliminated states to the highest 



eigenvalue for the retained states is large, a good approximation to the full-state linear 
model can be achieved. A property of the normal reduction technique is that the condi- 
tion X n = 0 forces the final value of the reduced models to the same steady-state value. 
s2 

This type of reduction is known as a type II method. An alternative method whereby 
Xg = 0 is identified as a type I reduction technique. 

The real-time, hybrid computer simulation of the Pratt & Whitney FlOO-PW- 100(3) 
turbofan engine was used to determine the linear models at the sea-level-static condi- 
tion. The state, control, and output matrices of equations (1) and (2) were generated at 
several power conditions along the operating line. The system matrices for the linear 
model representation are presented in appendix C. 

The state and control disturbances were limited to insure that engine performance 
would be maintained over its linear operating range. Because of the low-level disturb- 
ance signals, the output variations were so small that they caused large errors in the 
sampling process. The errors are mainly due to resolution and inaccuracies of the re- 
cording devices. Studies can be conducted to establish the optimum magnitudes of the 
disturbance signals to obtain good correlation between nonlinear and linear models. 
However, because of the tremendous effort expended in generating system matrices, 
the magnitudes of the distvirbances were maintained at a constant value. The model- 
reduction technique, described previously, was used to reduce the full-state (16th order) 
linear’ models to lower-order models. The system matrices for several reduced-order 
models at the sea-level-static intermediate power condition are given in appendix C. 

DISCUSSION OF RESULTS 
Nonlinear and Full-State Models 

Figures 6 and 7 present the comparison of nonlinear and full-state-linear-model 
responses to step changes in exhaust -nozzle area and main-burner -fuel flow. The non- 
linear responses were obtained from the F100-PW-100(3) real-time hybrid computer 
simulation. 

Figure 6(a) presents the chaise in fan speed Nj^ to a step change in exhaust-nozzle 
area. The two responses are at first quite close, but a slightly lower final steady-state 
value is obtained for the nonlinear model response. This error is about 12 rpm or 
11 percent. Similar responses for compressor speed Njj are shown in figure 6(b). The 
linear and nonlinear models exhibit different performance characteristics in that the 
final value of the linear model is much smaller than the nonlinear. The dynamic differ- 
ence between the nonlinear and linear models is caused by low signal levels. The ac- 
curacy and resolution of the analog -digital converters may be insufficient to detect the 
low signal levels. The results for fan turbine-inlet temperature T41 are given in 
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figure 6(c). The nonlinear and linear models exhibit similar characteristics. The non- 
linear model experiences a higher temperature overshoot (about 1 K) than the full-state 
model. The final value of the nonlinear response is slightly lower (about 2. 2 K). The 
extremely small magnitude of the final value (0. 15 mV) for the full-state response is 
causing inaccurate results. The responses illustrating the chaise in exhaust-nozzle 
pressure for a charge in exhaust-nozzle area are shown in f^ure 6(d). The char- 
acteristics of the two models follow closely with a final steady-state error of about 
0. 028 N/cm^ or 6. 6 percent. Thrust responses are shown in figure 6(e). Thrust is an 
output quantity so that the feed-forward term Du contributes to the initial thrust. The 
feed-forward matrix D represents the direct, linear relation between output y and 
input u. Slightly more undershoot can be noted for the nonlinear response. The final 
steady- state error is about 0. 17 kN. 

Figure 7 presents selected step responses to fuel-flow disturbances in the main 
burner. The change in fan speed Nj^ is given in figure 7(a). Althoi^h the nonlinear 
and linear models exhibit similar dynamic characteristics, a steady-state error of 
11 rpm, or 15. 9 percent, was obtained. The responses for compressor speed Njj are 
compared in figure 7(b). The two models exhibit similar characteristics, although a 
final steady-state error of 15 rpm, or 30 percent, was observed. The magnitude of the 
signal levels for the full-state model is small (about 2. 4 mV max) causing large errors. 
The accuracy of the analc^ digital converters is evident for the full state response. The 
results for main burner pressure for a fuel-flow disturbance are presented in fig- 
ure 7(c). The nonlinear response exhibits the same trend as the linear response; how-, 

o 

ever, a final steady-state error of 0. 5 N/cm or 17 percent was noted. The responses 
for fan turbine-inlet temperature T^j^ are shown in figure 7(d). The nonlinear response 
shows slightly more temperature overshoot and a 0. 8 K or 10 percent higher steady- 
state value than the linear response. The chaise in exhaust-nozzle pressure to a fuel- 
flow disturbance is presented in figure 7(e). Although the nonlinear response follows 
the same trend as the linear model; it had a final steady-state error of 0. 035 N/cm , 
or 14 percent. The change in thrust is given in figure 7(f). A higher final steady-state 
value (about 0. 13 kN or 15 percent) was obtained for the nonlinear response. 


Full-State and Reduced-Order Models 

The responses for full-state and reduced-order linear models for step changes in 
exhaust-nozzle area and main-burner-fuel flow are presented in figures 8 and 9. Step 
responses for several reduced-order models were generated to establish the model 
order that would most closely match the full-state model. The adequacy of the full- 
state and reduced-order linear models to describe the nonlinear engine was determined 
by applyii^ disturbances (area and fuel flow) to the various linear models and observing 
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the transient behavior. 

The change in fan speed to a step chaise in exhaust-nozzle area for several 
linear models is presented in figure 8(a). Excellent agreement was obtained for the 
full- state and seventh- order linear models. A slightly faster rise time can be observed 
for the lower order models (fifth and third). The dynamics of the linear models are 
nearly identical, and fan speed could be approximated by a lower order model. Similar 
step responses for compressor speed are shown in figure 8(b). Excellent agreement can 
be observed between the full- state and seventh-order models. A slightly faster rise 
time was noted for the fifth- and third-order models. The response characteristics of 
the linear nuodels differ from the nonlinear model (fig. 6(b)). As mentioned previously, 
the results for the linear models may contain inaccuracies due to the signal levels being 
too low for accurate detection by the conversion devices. Again, the full-state model can 
be approximated by a lower order model. The results for burner pressure (fig. 8(c)) 
show extremely small burner-pressure variations for a step change in exhaust-nozzle 
area. Excellent agreement can be observed for the full-state and seventh-order models 
and for the fifth- and third-order models. A slightly larger pressure undershoot exists 
for the two lower order models. The responses for fan -turbine -inlet temperature T 4 I 
are given in (f^. 8(d)). This state is a reconstructed quantity for the third-order model, 
as is evident from the initial temperature rise. The feed-forward term in the state 
output equation (29) contributes to the initial temperature rise. The magnitude of the 
signal levels (about 0. 5 mV) is extremely small and can cause difficulties in the samp- 
ling process. The responses illustrating the change in exhaust-nozzle pressure P.^ for 
an exhaust nozzle disturbance (fig, 8(e)) show excellent agreement between full-state and 
seventh- order linear models. Since exhaust-nozzle pressure is a reconstructed state 
for the third-order model, the feed-forward term H causes an initial pressure drop 
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of 0. 59 N/cm . After the initial transients, good agreement is obtained between the 
fifth- and third-order models. The thrust Fj^ responses are presented in figure 8(f). 
Because thrust is an output quantity, the feed-forward term Du causes an initial change. 
The agreement for the full-state and seventh-order model is favorable, but slightly more 
imdershoot is obtained for the full-state model. The third-order model response ini- 
tially differs from the fifth-order response, but after 25 msec follows closely. 

Figures 9 present selected responses to main -burner -fuel-flow disturbances. The 
change in fan speed Nj^ results (fig. 9(a)) show excellent agreement between the vari- 
ous linear models. However, the seventh-order model has a slightly higher rise time 
than the other reduced-order linear models. Since the dynamic characteristics of the 
full-state and reduced-order models are extremely close, the full-state model can be 
approximated by a lower order model. Excellent agreement can also be noted in the 
responses for compressor speed Njj (fig. 9(b)) for the full-state and seventh-order 
models and the two lower order models. The results for burner pressure P^ (fig. 9(c)) 
show that the responses for the various linear models are extremely close and allow the 
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full-state model to be approximated by a lower order model. Fan-turbine- inlet tem- 
perature responses are shown in figure 9(d). The temperature overshoot of the 

reduced order models is slightly higher than the full state model. However, good agree- 
ment between the various linear models can be observed. The fan-turbine-inlet tem- 
perature is a reconstructed state for the third- order model as indicated by the initial 
temperature rise. The change in exhaust-nozzle pressure P,^ for a fuel-flow disturb- 
ance is given in figure 9(e). Good agreement can be noted for the full-state and 
seventh-order models. A slightly faster rise time is exhibited by the two lower order 
models. Since exhaust-nozzle pressure is a reconstructed state for the third-order 
model, there is a small, initial pressure decrease in the response. The change in 
thrust is given in figure 9(f). Good agreement can be noted for the full-state and 
seventh-order models; whereas a slightly faster initial rise time was observed for the 
fifth- and third-order models. 

The eigenvalues for several full-state linear models along the sea-level-static, 
normal operating line are presented in table I. The eigenvalues for the full- state - 
linear and selected reduced- order -linear models are tabulated in table II. A graphical 
representation of the eigenvalues for the full-state and fifth-order models is given in 
figure 10. This analysis indicates that the important eigenvalues are located at (-4. 86, 
-6. 49) and (-4. 25, -6. 85). Similar results were obtained for the seventh- and third- 
order reduced models. Without further investigation these results imply that the dy- 
namic characteristics of the full-state model can adequately be represented by a second- 
order model, such as the two rotor speeds. However, additional system states may be 
required for engine control and protection. 


SUMMARY OF RESULTS 

A real-time, hybrid-engine simulation was formulated to evaluate the multivariable- 
control concept over specified power and flight conditions. The simulation was also 
used to generate linear models that describe the dynamic performance for low-level 
control disturbances. The latter task was performed to demonstrate linear-model- 
generation techniques and develop computer methods for sampling and recording the 
disturbance data. The linear -model matrices were computed on the EBM 360 computer. 

The linear models, valid about a nominal operating point, were used to obtain the 
response characteristics for exhaust-nozzle area and main-burner-fuel-flow disturb- 
ances. The linear models were cross-plotted with the nonlinear responses. Differences 
in the dynamic responses and final steady-state values can be observed for the linear and 
nonlinear models. The nonlinear models exhibit, for several conditions, slightly larger 
overshoot or undershoot. The difficulties encountered in the sampling and recording 
process are evident by the final value errors. In many cases the signal levels are ex- 
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tremely small and below the resolution of the analog -digital conversion modules. 

For large systems, such as the FlOO turbofan engine, the order of the state vector 
is large and may contain high-frequency information not required for control purposes. 
Thus, reduction of the full-state model to a reduced, lower order model could be a vi- 
able approach to simplifying the control design procedure. A reduction technique was 
developed to generate reduced-order models. 

The comparison indicates good agreement for full-state and seventh-order models. 
The initial dynamic characteristics of the lower order models differs from the higher 
order responses. The error can be attributed to the type of reduction and could pos- 
sibly be minimized by proper state selection for the reduced-order models. The type n 
reduction method, due to the Xg = 0 approximation, forces the final steady-state 
values to identical magnitudes. 

The advantages of the hybrid simulation are the ease and flexibility in setting flight 
and power conditions of the turbofan engine. However, the accuracy and repeatability 
accomplished by a digital simulation cannot be matched with the hybrid simulation. As 
noted by the nonlinear and linear results, errors can be introduced by the low signal and 
high noise levels. These small signals can also cause errors in the sampling and 
recording processes. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 16, 1978, 

505-05. 
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APPENDIX A 
SYMBOLS 


A state matrix n x n 

A transformed state matrix n X n 

An permuted state matrix n x n 

Ag reduced-order state matrix m x m 

Ajj partitioned state matrix m x m 

A ^2 partitioned state matrix m x (n - m) 

A2 j partitioned state matrix (n - m) x m 

A22 partitioned state matrix (n - m) X (n - m) 

B control matrix n x p 

B transformed control matrix n x p 

B^ permuted control matrix n x p 

B_ reduced-order control matrix m x p 

Bj partitioned control matrix m x p 

B2 partitioned control matrix (n - m) x p 

C output matrix r x n 

C transformed output matrix r x n 

permuted output matrix r x n 

C reduced-order output matrix r x m 

Cj partitioned output matrix rxm 

C2 partitioned output matrix r x (n - m) 

D feed forward matrix r x p 

D reduced-order feed forward matrix r x p 
s 

E state transformation matrix n x n 
Ep state permutation matrix n x n 
H state output matrix n x n 

permuted state output matrix n x n 
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Hg reduced-order state output matrix n X m 

Hy reduced-order state feed forward matrix n x p 

Hj partitioned state output matrix n x m 

Hg partitioned state output matrix n x (n - m) 

2 

P pressure, N/cm 

gas constant of air, 2. 87x10^ N • cm/^ K 
T temperature, K 

u control vector p x 1 

3 

V volvune, cm 

W stored mass, kg 

X state vector n x 1 

Xsj state vector (slowly decaying states) m x 1 
Xg2 state vector (rapidly decaying states) (n - m) x 1 
Xj transformed state vector n x 1 
Xg permuted state vector n x 1 
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appendix b 


HYBRID COMPUTER LISTINGS 


State Variable Initialization Subroutine (State 1) 


SUBkOU'l S lH iE::! <1 > 

SChI kr> BkHB TlON VHk J < Ir' i I ■' 

D I Mb- r-4S 3 ON HRDk C > 

COMMON/STH I kl/VHk3 .. !:■ 

DAT R RRDR < 1 > .. RRl Jk 1-' > .• HRi >k '■! !< > .• RRI 'k f' 4 • HHl ?k (. 1. > > .. RRI .'k < R .. RRI :*k V > 
± fiflDk <! & RRlJR < 9 > .. HRI :'k i.. :1 > • RRi ik (. J 1 . RhI :>K (! :) ?' > .. RRf ik t; -1 '< > .. RRl'iiR < J 4 

2 RRDk < 1 5 > .. RRDk t i R > .. RRI >k J ? > HRI -'k < :'l R > .. RRI :<k < j 9 > .■ RHl Jk 2R •• .■ RRi Jk C 

3 flflDk'l.22.Jj RRlJk':!2!.<>.-'4HRi.'Rki.. 4HHt-iRl.i,. 4HRRHl.i.. 4HR.1 J !.i,. 4HR:1 Rfi,. 4HR(-i8bi 

4 4HR0fc!5.' 4HR.1 -1 0.- 4HRti-'_ill',. 4HHij’i;'l fci.. 4HHi.iR?,. 4HRR.'-'('i.i 4HRi'ifeivi,. 4HRCHvi7.. 

5 4HR01?^ 4HR0?:0.. 4HR001 4HRi;:i:1 l.i.. 4HR0R0,. 4 Hh 0;>R.. ‘IHRtilil,... 

6 4HR105.. 4HR011 
C*****PLRCE Rr-4RI.DG;-. ] 1-4 Hul !.• 

CflLI.. t.!SC<2,. Jbkk> 

CALL 0SH< llr:kR.-:i 

C***:+:*REAU (.;0b4&01..k ijr-4rr RHki, .ibjhp-. 

CALL L!9C<0.. Ibkk> 

CALL QSC<J.. IbkkJ 
DO 4 21 

CALL QRRSCARtJk':; 3 :j. VHR3 < 3 j . .1 bkk :• 

4 CONT3:NUb 

c+*:+=**reai:j cor-4sai..b iwo rmri jt-it-k"-. 

CALI. (-.iSCtO.i it.kklJ 
CALL QSC<2, Jfckk.'J 
DO r. 1=1. 12 

CALI- Qk AS < ARI.JR < ] > . VRR I < 3 > .. .1 b;kk > 

5 COt4T.INUE 

VAR 1 < 1? .'J =- VAR I < 3 > 

VAR I< n > =- VRR I<; 5 > 

VAR I < 10 3- =- VAR I c: 10 > 

VRR I C 12 ;j =- V AR I c; 1.2 > 

CALL QRAS < AADR ( 22 > .. VAR 3 < 22 > . J I- kk ;j 
CALL. QRAS < RADR C 23 > . VRR I. C 23 > . 1 Fkk > 

VAR I. < 22 ;j =-- VAR I < 22 > 

VAR I < 23 > =-- VAR 1 1 23 > 

C***++SET 1C POTS FOR ALL. INTEGRATORS 
CALL QWBDCS<42. VARKi;*.. IG. lERkiJ 
CALL QWDCS<58. VARK:221. 1. lERkU 
CALL QWDCS<71.. VARK:23 :j.. 1.. IFRRj 
i::* 4:4::+:+PLACE ANALOGS IN IC 
CALL QSCCl. lERRO 
CALL QSICdERRl 
C-^**^:* INITIALIZE FOR STATE 2 
L=0 

D=. 02S 

CALL QSCC0. IERR> 

CALL QSCCl. 1ERR> 

RETURN 

END 
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state Variable Disturbance Subroutine (State 2) 


hOk'IRHN CLiHHUfek kRV. I b V. Jfei© 


SUBki: iij I 3. 1-JE B 'l H Tt- 2 (. L > 

scBi ki;' FkHuiioN vfiki l:-.. vHk';:2;;0 
CGMMUH.-'S I BTIr 1 /VBk I .. D 
CFlLl Q!r*C »•! 1 F RR } 

C R L 1. Q S C < P .• J. F R R! !•* 

^ L=L+1 

IF<L . o r. 4:<:' GO I U 4i' 

IFCL. G T. > GO I II PGG0 



GO lO 

•:;i. 2. 7. 

4. 8. V'. 8. 

9. :1 0. :1 1 . 1 2. 

0ti0 

LL==:L-P 

1 




GO ro 

<22. 27. 

24. PO. 28. P 

7. 28. 29. 70. 


i 41. 42 

. 40 >. 1. 1.. 



00 0 

FORMRl 

<7K. 3 2> 



.+ + .+. 

+ NH LTIF 

•( UkPRrMi:. 

T 


1 

VRR 1 

-VRk 1 < J 

> •+*j.:‘:+.VRk J 

‘[i 


CRLL Q 

'WDl;s<4P 

'. VhR‘;.:3 X. X. 

ibRR:;- 


GO 10 

L • 0 



P 

VRR‘..;i > 

■ VHk.l *21 

■j — i;» +-'v’Hk .1 

> 


CRLl, (.! 

||AlDi;.:S<4P 

’. VHk*:. J ..1 . 

J b kk 


GO lU -.'O 

I. ;•+ + +■•+ +:NL 1 S I URRHI'-lL^h 
S vRR‘..;J > -vRRx ':..1 > 

VRRC^tD-VRhO 'U-'^+b + VHk 1 'vc’> 

CRLL (.!WBDC:r.<4P. VRK*::! P. J.bR:R> 
GO I C.» Lt0 

4 VRR P > VRR J < P > ’“O += VHR i <y> 

CRLL QI'-IDCSC. 4I<.* VRR<P“>.> !.■ I FRk> 
GO i U 00 

1 ..+..+..+::+::+:W41 1.:* I STURRHh4C:F 

5 vrr<p>^-vrri<p:;‘ 

VRR VRR U. 'S > +D=+ VHk i < :< > 

CRLL 0WBDCS<43^ VRR':;P>.. P/ ibRR> 
GO ro 0© 

6 VRR :< > ^VRR 3 3: > -D+VHR i > 

CRLL i;:-!WDC:S*:;44.. VRR<’<>;* I.* 1FRR> 
GO TO 00 

C + D I SI URBRr-4Ct: 

7 VRR VRR 1<2> 

VRR < 4 > —VRR I < 4 > +D* VRR J < 4 > 

CRLL QWBDCS<44.» VRRCjP^.* P^ 1ERR> 
GO "10 00 

8 VRR c: 4 > ^VRR i < 4 > ~I.:«+VRR 1 C 4 > 

CRLL QWDCS<40. VRRC4::'.. J .. ltRR> 
GO TO 00 

C + + + ++P41 i::* I STURBRf-4CE 

9 VRR 4 > - VRR i< 4 > 

VRR C ^ VRR I O +D+- VRR J 0 > 

CRLL. QWBDCS < 40.. VRR c! 4 > .. P:.. 1 F.kR > 
GO TO 00 

10 VRR<0>^VRRI '::o::‘-L>^^:VRRI<0> 

CRLL 0WDCS»::46. Vf\RCJ>.> 1.- IlrRR> 
GO TO 00 


J 4. O. 6. 1 7, J 8.. 9. 20. 2: 

!:C2. 74. 70. 78. 77. 78. 7: 


U, L 
S'. 40. 
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inmiiiiiiiiiiiiiiiiiiiiii 


CH- *++:+:W6 if 1 STUkBflr-4Cf- 

VflR ■; 6 y - VRR I < 6 > +D*VRR 1 1 6 > 

CALL QWBDCS<46.. VRR<;5>, 2, IERR> 
GO TO 00 

12 VRR < 6 - VRR I < 6 > -D*VRR I < 6 > 

CALL QWDCS<47, VRR<6>^ ±. lERR') 
GO TO 00 

L;:+.+>t:**WG6 D 1 STURBRr4CE 
1^-: VRRc:6>-VRRI<6> 

VRR C ? > -VRR I C 7 > +D*VRR I < 7 > 

CALL QWBDCS<47.. VRR<6:), 2, IERR> 
GO TO 00 

14 VRR •; 7 > - VRR jU 7 > -L):+: VRR 1 1 7 > 

CALL QWDCS<:48^ VRR<7>^ 1.. I ERR) 
GO TO 50 

C=+.++:*>t:P6 1 .) I ST URBRr-4CE 

15 VflR 7 ) =- VRR I < 7 > 

VRR < S ) -VRR I < B > +D+= VRR I < 8 > 

CALL 0WBDi;.;S<48.. VRR<:7)^ 2. I ERR) 
GO TO 00 

16 VRR < 8 ) —VRR I < 8 ) VRR KB) 

CALL QWDCS<;49, VRR<S).. Ij I ERR) 
GO TO 00 

i:;;:+::+:+*h<W7 D I ST URBRrTCE 

17 VRR 8 ) -VRR K 8 ) 

VRR c; 9 ) - VRR I < 9 ) +Dh-:VRR I < 9 ) 

CALL QWBDCS<49j. VRRCS).. 2.. lERR) 
GO 10 00 

18 VRR 9 ) -VRR K 9 ) -D*VRR 1 C 9 ) 

CRLl., QWDCSCOO. VRR<9)^ 1^ I ERR) 
GO TO 00 

DIST URBRfTCE 

19 VRR < 9 ) - VRR I c; 9 ) 

VRR <; 10 ) =VRR I 10 ) +D+VRR K J 0 ) 
CRLl. QWEDCS C 00.. VRR < 9 ) .. 2.. lERR) 
OO T O Oi0 

20 VRR < 1 0 ) = VRR K 1 0 ) -D* VRR K 1 0 ) 
CRLL. £!WDCS';;OJ... VRR (.10).. 1.. I ERR) 
GO TO 00 

il:*+*:+:*T12 !.:■] STURBRfTCE 

21 VhR C 10 ) VRR 1 1 10 ) 

VRR 11 ) =.VRR 1 <. 1 1 ) +t)4:VRR 1 < 11 ) 
CRLl. QWBDCSCOl/ VRR<10).. 2.. lERR) 
GO TO 00 

22 VRR < 11 ) =VRR 1 < 11. ) -D+VRR KID 
CRLL 0T4DCS<02.. VRRCll ).. 1 .. I ERR) 
GO TO 00 

C:+ D I STURBRNCE. 

2s VRR < 1, 1 ) -vm. K 1 1 ) 

VRR •; 12 ) =VRR 1 < 1 2 ) +D*VRR ] C J 2 ) 

chll i;-.!Wedce;(:52.. vrrcii ), 2 .. jerr) 

GO TO 00 

24 VRR c 12 ) =VRR K 12 ) -D* VRR JC 1 2' ) 
CRLl.. UWDCSCOS.. VRR<12).. 1. lERR) 
OO TO fO0 

C++*+=+l.43 0 1 ST URBRNCE 
2.0 VRR C 12 ) = VRR I c; 12. ) 

VRR C 13 ) -VRR K 1 3 ) +D*VRR K 1 3 ) 
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CRLI. c**.- Ifr*RR> 

GO I O 00 

26- VRR :1 > -VRR 1 < 1 2 > RR 3 < :1 !:> > 

CRLI . C!WL:'*1;S < ‘-*4.. VRR < ^ ::=• > -• J -• i R RR > 
GO I Cj 00 

D1STURBRNC6 

27 VRR <! :11:> > - VRR J •• ^ > 

VRR '■■: 14 > = VRR J < 1 4 > +D'+*VRR 1 < 1 4 > 
CRLL. QWBDl:R < 04.. VRR < i > .* 2. 1 B RR > 
GO 10 00 

28 VRR < 14 > =VRR 1 1 4 > VRR I < 14 > 

CRLL QWDC8 < 00.^ VRR < J 4 1 ^ I ERR > 

GO rO 00 

i;. 5+:p 4 l.>3. S ( URBRr‘4CE 

29 VRR < 14 > —VRR! I < 14 > 

VRR 'v Is.' > - VRR I < J 0 > +D+VRR 1 < l.Oi > 
CRLI OWBDCS c; 00.. VRR < 1 4 > .. 2.. 1 R RR > 
GO 10 00 

30 VRR < lOt - VRR I C 1 0 > -L:*=+- VRR 1 < 1 0 > 
CRLI. OWDCS < 06.. VRR i1 0* > .. J .. 3 R RR > 


GO I O 00 

t;:+..+ >K.:+: + |..44 ILl S 1 URBRNL-It 
31 VRR !*-• - VRR 3 J 0 > 

VRR 16 — VRR I J 6 > +L':+-VRR 3 i1 6 .> 

CRLL. OWBDi; :R 06.. VRR < i > .• 2.. 1 L RR > 

GO f L« !-«0 

22 VRR :1 6 > - VRR 3 J 6 > -L>:+*VRR: 1 *• -1 ^ 

CRLI uRDi:;S 0?.. VRR 1 6 > .. -1 •• 3. L RR > 

GO TLi 00 

(.;.+• +.:+;-+:+l--4F 4 L> .3 6 3 lJRBRr‘4C.:L 
2 R R J 6 > - V R R 1 < 1 6 

CRLL UWDCS c: Li 7.. VRR < J 6 > -• 3 -• I L RR > 

C R L. I L! 6 C < 0 .• 3 L. R R > 

C R L I . y. 6 C < :1 1 1. R R > 

VRR < :1 7 .i' -D + VHR I < 1 7 > 

CALL. QWJOH;-.<:VHk<d.V’>.. 13.. 

GO I LI L'O* 

...4 wF4t‘1 , k-i3.SV'iJS.-‘';. > 

IFCVHkl':;!?::'. t I.Wk4M> iiO 10 lOfJW 
C R L I y .6 C c‘ 0 .• i L R R > 

CRLL I’-! 6 C: •:! ;1 .. 3 1 h! h! > 

R R 1. 7 > - - D •+• 'v' R R 3. J V ' > 

CRLI. L'W JDRS < VRR c: J 7 > .* -t 2.. J Ir RR > 

GO I L-I OOi 

|“ :+:4: +::+ :+4.4F 7 D 3 6 T URBRf ^CL 

25 CRLL. C!SC':;0.. 1ERR> 

CRLL WBC':; J .* 3 ERR;:« 

VRR<17>=. 06 
VRR'::18>=L' 

CRLL QW JDRS < VRR < 1 7 > .. l.l.* 1 ERR > 
CRLL QW JDRS < VRR < IS > . ^ 4. IE RR 
GO TO 05 

26 IF<VRRI'::i8>. L.T. D> GO TO 1000 
CRLL QSC<0.* lERRC' 

CRLL QSC<1.. IERR> 

VRR<1S>='~D 

CRLL QWJDRS<VfiR<lS:L. 14. IERR> 


I INIVIIIIIIIIIillllllllllllllllll 


GO TO r35 

i;. X SI URBRT4CE 
3? CRLL. QSC<0.. 1ERR> 

CflLl 0SC<1. 1ERR> 

vflR<iw;j=. 0s 

VflR < 19 > --L)*VRR I < 19 > 

CRLL QW JDRS<VfiR<18>. 14^ ].ERR> 

CALL OWJDRS<VRR<19>,. 3 5, 1HRR> 

GO TO 55 

38 CRLL QSCCfei, 1ERR> 

CRLl OSCCl, 1ERR> 

VflR < 19 ;■ -D+ VRR I C 19 > 

CALL. OWJDflS < VRR :1 9 > , i 5, ] ERR > 

GO 10 55 

C:+:if;*4::+:C I V V I :• ]. STURfciRrTCE 

39 CALL C-!SC-:;0, lt.RR> 

CflLl OSCCl, 1ERR> 

VflR<;i9>-. OS 

CflLl OWJDRS < VRR C :1 9 > , :1 Ci, 1 E RR > 
IF‘:;VRR].<i:-'0>. L1 , 99999S+lj> GO '| O :1 080 

VflR <■ RO > =— L:> 

CflLl Ot'J.J DRS VRR < 20 > , 3 8, X E RR > 

GO 10 !.'5 

40 IF<VRE:l <R0>. GT. — L5> Oil TO 3 000 
CflLl OSC<0, IERR> 

CflLl OSCCT, XERR> 

VflR <! R(-i 3 =D 

CFlL.i . OWJDRS < VRR C 20 > , 3 8, J E RR ;■ 

GO 1 O 55 

i;,*.+.:+::+:+:RCVV GIST LIRBRrTCb 

41 CRLl. OSC'iO, lERR.:- 
CflLl. l.•.!SCtl, ]ERR> 

VflR';.::-:Ei>=, os 

CflLl. OW.J DRS C VRR < 2fei > ,3 8, ] E RR > 

IF< VRRI <:23 >. I I , 99999S+Ij> GO Ri 3 000 

VflRc.;:;;1.>=-G 

CflLl. 0W.J15R.S < VRR 23 > , 3 ?, 1 E RR :> 

GO I O 55 

42 IFt. VRRl <23 >. G'f. -IT;' GiO f ti 3 000 
CflLl 0SC<0, JERR> 

CALL 0SC<3 , 3ERR.;i 
VflR <23 >=D 

CflLl OW.J[.:'R,S < VRR < 2:1 > , 3 1 E RR .J 

GO I (.1 55 

45 CflLl. 0SC<0, iERR> 

CALL OSC'rll., 1ERR> 

VflR < 2'1 3 =. OS 

CflLl. OW.JDRS < VRR < 23 3 , 3 V.. I ERR 3 
CflLl.. Q;SC<2, 1ERR3 
CflLl , OSOR< I ERR 3 
4? CflLl 0SC<0, 3 ERR 3 
CALL USC<3, 1E.RR3 
GO I O 80 

50 CflLl. OSC<0, 3 ERR 3 
CflLl. USCCl, .IERR3 
35 T VPE. 3000, I . 

CflLl OWCIJ. <1, TRUE.. , 3ERR3 
60 RETURI-4 
END 
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APPENDIX C 


FlOO SERIES n-3 SYSTEM MATRICES AT VARIOUS POWER CONDITIONS 


Full-State System at Sea-Level-Static - Intermediate-Power Condition; 
Power Lever Angle, 83° 


MATRIX A 



1 

2 

3 

4 5 

6 

7 

8 

1 

-4.721 

0. 1223 

-0. 2263 

-8629. -1.588 

728.7 

-3274. 

-35.01 

2 

-0.1671 

-5.230 

191.7 

-22,44 0.0000 

-241.5 

-377.4 

-3.868 

3 

0.0000 

0.4278D-01 

-ice. 8 

0.0000 -1. 122 

0.4669 

0.0000 

0.0000 

4 

0.7906D-03 

0. 3357D-02 

0.1748 

-23.32 -0.9354D-02 

-0.6769 

0. 2422 

0. 3552D-02 

5 

-C.3012D-02 

-0. 2971D-02 

-0.6395D-01 

5.360 -250.8 

-0.7781D-01 

1 .93P 

-0. 2037>D-01 

6 

0.3464D-C1 

0. 1188D-02 

21.61 

1110. 0.7475 

-83.41 

7 . 740 

0. ^0e*7D-01 

7 

0. 1732D-01 

-0.2567D-01 

0.8527D-02 

-0.5360 O.OCOO 

0.3113D-01 

-1.35*7 

-0. 16240-01 

8 

0.21 08D-02 

-0. 1188D-02 

0.0000 

0.0000 0.0000 

0. 3113D-01 

888.3 

-10.51 

9 

-0.7518D-05 

0. 9030D-03 

-4.360 

0.1608 0.5607D-02 

0.1245D-01 

0.4847D-01 

-0. 101-70-03 

10 

0.6325D-03 

0.2163D-01 

4. 184 

-0.6432D-01 -0.4488D-02 

-0. 3113D-03 

1.244 

0. 14 1-70-01 

11 

0.1852 

0.4753D-02 

0. 12-79 

5.360 0.0000 

0.1556 

•722.9 

*7.835 

12 

0.7364D-01 

0.7486D-02 

0. 6395D-02 

-589.1 -0.5610D-01 

U9.88 



13 

-0.7590D-03 

0.3541D-03 

O.OOCO 

24.45 0.7480D-03 

0.6894 

0. 1550D-01 

0. 40600-03 

14 

0.6795D-05 

0. 1070D-C4 

-0. 3840D-03 

-0.2410D-01 -3.180 

0.1402D-02 

0.0000 

-0.91260-04 

15 

0.0000 

0.2675D-02 

O.OCOO 

12.05 -2172. 

0.0000 

-4.358 

0.0000 

16 

0.3765D-01 

0.2970D-02 

O.OOCO 

20.10 0.0000 

-0.7781 

0. iret’D 05 

115.8 



9 

10 


11 

12 

13 

14 

15 

16 

1 

91.75 

-14.73 


-0. 1 1C4 

-610.0 

40.65 

33.34 

0.2871 

0.1837 

2 

-0.2691D 05 

-0. 14220 

05 

-13.54 

1.501 

26.73 

0.0000 

0 . 1 P 1 7 

-0. 3627 

3 

6570. 

0.2937D 

05 

22.96 

-1.075 

0 .OOOC 

0.0000 

-0.202° 

O.OCOO 

4 

97. 97 

-2.430 


-0.20800-02 

0.53760-01 

0.0000 

0.5889 

-0. 1692D-02 

0. 0000 

5 

-14,41 

0.0000 


0.0000 

153.7 

-2732. 

1370. 

28.85 

0.861=90-01 

6 

-5821. 

-69.42 


0.3744 

1 . 433 

0.0000 

15.70 

-0.6*7540-01 

0.3464 

7 

C.OOOC 

16.67 


0,10400-01 

-C. 71670-01 

1. 277 

O.OCOO 

0.0000 

-0. 9871 

8 

17.29 

8.333 


0. 0000 

0.5017 

3.830 

-6. 281 

-0. 5411D-01 

-0. 1732 

c 

-101.9 

1468. 


1.078 

0.17950-02 

0.63840-01 

0.0000 

0.67640-03 

0.85450-03 

10 

-0.57630-01 

-1483. 


-1. 089 

0. 71670-02 

0. 38300-01 

0.1099 

-0. 67640-03 

0. 51950-02 

1 1 

0.0000 

7833. 


-14.27 

0.3584 

0.0000 

-7.852 

0.0000 

0.0000 

12 

4. 322 

-3.472 


0.21840-01 

-1C3.5 

2268. 

-1. 570 

-23.93 

12. 18 

13 

0.57630-01 

-0.2500 


-0. 14560- 03 

-0.48730-01 

0.51070-01 

0.7852D-01 

-1.006 

0.9971 

14 

-0.1296 

0.1249 


0.70220-04 

-0. 16110-02 

0.0000 

-66.90 

1 .00*7 

0.7*"86D-03 

15 

-129.6 

-31.23 


0. 23380-01 

2071 . 

14.36 

-159.0 

0.4566 

-0. 3894 

16 

-1045. 

173.6 


-0. *7 *7990-01 

-1954. 

-63.84 

39.26 

2.621 

-25.54 


MATRIX B 

1 

2 

3 

4 

5 

1 

-It .21 

1 .489 

O.OCOO 

0. 0000 

0.36910-02 

2 

6.539 

-262.9 

-58.01 

-0.30800-02 

0.0000 

3 

-1. 874 

0.0000 

-16.62 

0.2559 

-0. 13040-02 

4 

0.46840-01 

0.21C4 

0. 1386 

-0.36750-04 

-0.10870-04 

5 

0. 3123 

-0.3507 

-2454. 

0. 1471D-02 

n. 1135 

6 

0.0000 

6.312 

-5.541 

-0.14710-02 

-0. 8694.0-03 

V 

C. 9367 

-C. 9118 

-1. 108 

-0. 44120-03 

-0.43470-04 

e 

0.1249 

0.2806 

1. 108 

-0.14710-02 

-0.43470-03 

9 

-0.62450-02 

0. 45600-01 

0. 1386 

0.26470-03 

-0.26080-04 

10 

0.37470-02 

0.9427 

0. 1330 

0.11760-04 

0.12170-04 

1 1 

0.6245 

4.209 

0.0000 

0.0000 

-0.26080-02 

12 

-0.3435 

0.5962 

3.047 

-0. 22060-03 

-0. 17390-03 

13 

-0.49960-02 

0.23850-01 

0. 1 1 080-01 

0.86230-05 

-0.69560-05 

14 

-0.28080-02 

0.31540-02 

-84.06 

-0.32940-05 

C. 2*7390-03 

15 

O.COOO 

1.570 

-2144. 

-0.66170-02 

-0.29340-02 

16 

16.39 

11.40 

-8 3.11 

-0. 91910-02 

-0. 54320-03 
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MATRIX C 



1 

2 

3 

4 

5 

6 

8 

1 

0.1860D-01 

-0.9506D-04 

0. 0000 

-0.3431 

0. 1 19'?D-01 

-0.2489D-02 -0.c3C2D-C1 

0. 0000 

2 

0.0000 

0.7842D-05 

-0. 1876D-03 

-0. 1175D-01 

-C.8229D-C3 

0.0000 O.OOGC 

-0.8931I)-0ii 

3 

C.OOOO 

-0.5703D-02 

O.OCOO 

8. 576 

535. 6 

o.croo o.oocc 

c. c^cc 

4 

-0.2412D-03 

-0. 3801D-C3 

8. 165 

-1.717 

O.OCOO 

0.0000 0.62C0 

-0. 3251D-02 

c 

-0.1517D-06 

0.6074D-07 

-0.5322D-05 

-0 .26630-03 

-0. 2811D-04 

-C.39"8D-05 -0.20t*0D-01 

-C. 11ii2D-03 

6 

0.1337D-06 

0.0000 

-0. 122GD-07 

0.0000 

O.OCOO 

O.OCOO -0.8ii28D-0U 

-0. 91 18D-0«= 

7 

0.5745D-09 

0.4079D-08 

0.0000 

0.0000 

0.0000 

0.0000 0. IICcd-05 

0. 1162D-07 



9 

10 

1 1 

12 

13 

m 

15 

16 

1 

O.OCOO 

0.4445 

0.3328D-03 

0. 1 147D-01 

C. 1021 

O.OCOO 

0. 1082D-02 

0. 0000 

2 

0.0000 

0.0000 

0.0000 

0.0000 

O.OCOO 

O.OCOO 

-0. 1i’P8D-03 

0. OOCO 

3 

-138.3 

44.44 

C.4S92D-01 

-1. 147 

0.0000 

55 2.8 

C.OOCO 

0.0000 

4 

-5560. 

2.224 

O.OCOO 

-0.5726D-C1 

O.OCOO 

1.258 

O.OCOO 

0 . 

5 

-0.7365D-03 

0.3550D-03 

-0. 2499D-06 

-0. 26S7D-C4 

-0. 1632D-03 

0. = 7810-03 

-C. 5CP3D-05 

0. 1-762D-02 

6 

0.0000 

0.0000 

o.oooc 

C.OOOO 

O.OCCO 

O.OCOO 

O.OCOO 

0.0000 

7 

C.OOOO 

-0. 2252D-04 

-0. 1 686D-07 

O.OOCC 

O.OCOO 

C.OOOO 

0.0000 

O.OCCO 


MATRIX D 

1 

1 0.6894 

2 

-0.673aD-01 

3 

0.8860D-01 

4 

O.OCOO 

5 

0.41"3D-04 

2 

0.0000 

1.000 

O.COCO 

C.OOOO 

-C. 28690-05 

3 

1.998 

-2.245 

5754. 

2353D-02 

0. 69560-03 

4 

0.0000 

0,0000 

-0.8852 

-C.2349D-03 

0.0000 

5 

0.0000 

0,34*’9D-04 

0. 1416D-03 

0.0000 

-0. 21560-07 

6 

0.0000 

0.0000 

O.OOOC 

0.0000 

0.0000 

7 

0.0000 

0.0000 

0 . cooo 

O.OOOC 

O.OCOO 


Full-State System at Sea -Level -Static Condition; Power Lever Angle, 70® 


MATRIX A 



1 

2 

3 

4 5 

6 7 

8 

1 

-3.754 

-0. 27800-01 

-0. 1726 

-7600. 0.0000 

713.0 -2386. 

-24. 35 

2 

-0.1495 

-3.806 

200.7 

36.94 -1. 890 

-23*7,0 -520.0 

-5.218 

3 

-0.48660-02 

0.18710-01 

-101.0 

0.0000 O.OCOO 

1.163 O.OCOO 

0.0000 

4 

0.72990-03 

0.25260-02 

0.1829 

-21.02 -0. 3385D-C1 

-0.7170 C.9P0? 

C. 87440-^^2 

5 

-0.81110-03 

0. 18710-02 

O.OCCO 

-5.880 -230.8 

0.2907 2.1"S 

O.OCOO 

6 

0.2'’580-01 

0.2370D-01 

19.13 

864.3 0.0000 

-75.36 17.42 

0. 1749 

7 

0. 15570-01 

-0. 13640-01 

0.0000 

-0.5880 O.OCOO 

0.19380-01 -6.753 

-0 . 74 320-0'’ 

8 

0. €4880-03 

-0.499CD-03 

-0. 2168D-01 

1,176 0.9025D-01 

-0.7”52D-01 693.8 

-12. 94 

9 

-0.32440- 04 

0.51150-03 

-4.463 

-0.2352 -C.1805D-01 

0.14^40-01 -0.54tt’D-01 

-P. =4670-03 

10 

0. 11940-02 

0. 16070-01 

4. 300 

0.9407D-01 -0.4513D-02 

0.38760-02 4.753 

0. 4 06 ID- '"I 

11 

0. 48660-01 

0.0000 

0.5420D-01 

O.OCOO O.OCOO 

-0.1938 222.2 

2.2*73 

12 

0.65530-01 

0.5C52D-02 

0. 1355D-01 

-453.9 0.9025D-C1 

45.50 1.525 

1.423 

13 

-0. 56770-03 

0. 27190-03 

0.5421D-03 

21.74 -0.902=0-03 

0.*7364 0. 82780-01 

0. 6?’^ = D-C3 

14 

0.0000 

0.28060- 05 

0. 1220D-03 

0.0000 -3.C5C 

0.13CPO-02 O.OOOC 

_n , 0 01 PD-C4 

15 

0.72970-02 

-0.28050- 02 

C. 1220 

13.23 -2138. 

0.0000 -4.9C1 

- 0 . 4 c 19D-01 

16 

0.0000 

0. 15590-02 

-0.6774D-01 

0.0000 -0.5646 

0,2422 9841. 

91. 29 


9 

10 

11 

12 13 

14 1 = 

16 

1 

-101.7 

102.5 

0.0000 

-588.3 7.183 

61,74 -0.3966 

3,060 

2 

-0.2809D 05 

-0.15370 05 

-13.47 

-1.849 14.18 

-1*7.41 -0.4 = 96 

0.7549 

3 

4982. 

0.2439D 05 

18.27 

O.COOO O.OCOO 

24.93 -C.2243 

-C.54 07 

4 

94. 21 

-0.4C24 

-0. 1 138D-02 

0.4414d-01 -0,3386 

0.2077 0,18680-02 

0.4=040-02 

c 

-31.94 

-8.050 

0.0000 

147.0 -2268. 

1276. 25.1° 

-0 . =4 C*7 

6 

-4120. 

16.08 

0. 3529 

2.206 6.762 

O.OCOO 0 . 

1800 

-7 

-3.194 

3.220 

-0. 1 138D-02 

O.COOO 0.6771 

1.662 -0.7U76D-02 

- 0. 0*7 3 1 

8 

31.94 

3. 220 

-0. 22770-02 

0.1766 -2.70° 

-3. 324 0.70 7^1 )- 01 

-C. 2884 

9 

-99 .00 

1356. 

0.9431 

0. 1324D-01 -0. 57710-01 

0.2008 0.3744D-03 

0. 360ar-02 

10 

-0.7026 

-1360, 

-0. 9458 

-0.15010-01 0.1625 

0. 33240-01 -C. 448=0-0^ 

C. 18C2D-^2 

1 1 

0.0000 

8001. 

-14,33 

O.OOOC -6.771 

-8.310 C. *^4760-01 

0. CCOO 

12 

-7. 984 

0.8050 

0. 22200-01 

-100.1 1883. 

-1. 247 -2C.92 

1 . 55 

13 

-C. 1916 

0.1288 

C. 68310-04 

-0.49430-01 -0.6771D-C1 

-C.6649D-01 -O.occi; 

1 .018 

14 

-0.2156 

0.1087 

0.76850-04 

0.0000 0.45700-01 

- 54.94 1 . 00 c 

C. PI C8D-'^3 

15 

-71.81 

36.20 

0. 25620-01 

208"'. 0.0000 

-3*7,4 0 -C.168? 

-0.4052 

16 

39.96 

301.9 

O.COCO 

-1549. O.OCOO 

353.2 C.1P68 

-179.0 
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KATEIX B 



1 

2 

3 

4 

5 

1 

-77. 53 

-4.519 

33.35 

-0.2272D-01 

-0. 18400-02 

2 

-2.614 

-182.8 

32.92 

-0.1345D-01 

-0.18150-02 

3 

0.0000 

2.129 

31.43 

0.2698 

0. 13000-02 

4 

-0.1560D-01 

0. 1242 

0. 1310 

0. 10710-03 

0. 21670-04 

5 

0.3122 

0.5323 

-1771. 

0.10710-02 

0, 1296 

6 

-0.6244 

2.483 

10.47 

O.OOOC 

-0.43280-03 

7 

0.9989 

-0.8871 

0. 5239 

-0.21410-03 

0.43340-04 

8 

-0.1249 

-0.4968 

-2.096 

C. 42820-03 

0.0000 

9 

-0.6244D-02 

-0.5326D-C2 

-0. 1 572 

0. 31050-03 

-0.86680-05 

10 

0. 1124D-01 

0.7750 

0.4191D-01 

-0.12850-04 

0.13000-04 

11 

0.6244 

0.3549 

0. OOOU 

-0.21410-02 

-0.21670-02 

12 

-0.9366D-01 

0.2661 

-0. 5239 

0.21410-03 

0. 65010-04 

13 

-0.3746D-02 

0. 1206D-01 

0.0000 

0.25690-04 

-0.34670-05 

14 

0.0000 

-0. 3194D-02 

-68.38 

0.0000 

0.27300-03 

15 

0.0000 

-1.597 

-1733. 

-0.48140-02 

-0.19500-02 

ie 

0.0000 

-0.4440 

6.555 

0.0000 

-0.75850-02 


MATRIX C 

1 

2 

3 

4 

5 

6 

7 

8 

1 

0. 17880-01 

0.39930-04 

-0. 17350-02 

0.0000 

-0.72160-02 

0.0000 

-1.290 

-0. 12940-01 

2 

0.17910-04 

O.COOO 

0.0000 

0.0000 

0.10110-02 

0.43410-03 

0.48800-02 

O.OOCO 

3 

0.51910-02 

0.19960-02 

0. 2602 

9.407 

520. 6 

-0.6202 

-3.486 

-0.34980-01 

4 

0.25980-03 

0.0000 

9. 176 

0.0000 

-0.2166 

-0.62000-01 

-0. 3489 

-0. 10490-01 

5 

0. 20760-06 

-0.15970-06 

-0.26C0O-05 

0.28200-03 

0. 0000 

-0.15510-04 

-C, 34300-01 -0.18850-03 

6 

0.13550-06 

-0. 11900- 0<5 

0.51700-08 

0.0000 

-0.43040-07 

0.18480-07 

-0.88710-04 -0. 91310-06 

7 

0.30940-09 

0.41640-08 

-0. 20680-07 

0.0000 

O.OCOO 

0.0000 

0. 12470-05 

0. 18760-07 



9 10 

1 1 

12 

13 

14 

15 

16 

1 

0.0000 -0.7728 

0. 36440-03 

0.0000 

0.0000 

-0.2660 

0. 1 1960-02 

0. 86500-0? 

2 

0.0000 0.69550-01 

0.4918D-04 

-C. 1 9070-02 

-0.151-^0-01 

-0.18620-01 

0. 16750-03 

-0, 40370-03 

3 

0.0000 -77.28 

-0. 54640-01 

0.0000 

10.83 

O.OOOC 

0.0000 

0.2884 

4 

-5401. 0.0000 

0. 54650-02 

0. 0000 

1 . 08 2 

3.990 

-0.23920-01 

O.OOCO 

5 

-0,20440-02 -0.25800-03 

-0. 20030-05 

-0.49460-04 

0.54180-03 

-0,66500-03 

0. 35860-05 

0.28600-02 

6 

0.0000 0.0000 

C. 21-710- 08 

0.0000 

-0.64580-06 

0.0000 

0.0000 

-0. I-'IQO-O? 

7 

0.12180-04 -0.27640-04 

-0. 19540-07 

0.0000 

0.0000 

0.0000 

0.0000 

-0. 69750-07 


MATRIX 0 

1 

2 

3 

4 

5 

1 

0.9690 

0.17030-01 

0. 0000 

O.COOO 

-0. 13870-04 

2 

0.0000 

1.000 

0. 1 1730-01 

0.0000 

0.0000 

3 

C.S990 

-0.5678 

3738. 

0.0000 

0.0000 

4 

-0.1990 

0.0000 

0. 0000 

0. 34210-03 

0.13860-03 

5 

0.99920-04 

0.0000 

0. 50310-03 

-0. 68610-0*7 

-0.27730-07 

6 

0,59550-07 

0.67680-07 

-0.49960-06 

0.0000 

0.41330-10 

7 

0.0000 

0.13540-06 

0. 19980-05 

0.0000 

0.0000 




Full-State System at fcJea- 

-Level-Static Condition; Power Lever Angle, 52® 


MATRIX A 









1 

2 

3 

4 

5 

6 

-7 

8 

1 

-2.465 

-0,68800-01 

0. 1448 

-2726. 

0.5928 

585.9 

-2013. 

-18. 84 

2 

-0,1846 

-1.835 

192. 7 

27.88 

5. 84P 

-222.5 

-552.0 

-5.086 

3 

-0.10580-01 

0.23340-01 

-93.34 

0.0000 

3. 351 

1.493 

29.56 

C.OOOO 

4 

0. 17640-02 

0.90770-03 

C. 1928 

-16.47 

0.2793D-01 

-0.95*79 

C. 9854 

0. 8758D-02 

5 

0. 17640-02 

0.64830-03 

-0. 34120-01 

-3. 328 

-216. 4 

-0.1244 

0.0000 

-0.7007D-01 

6 

0.58190-01 

-0.77780-02 

18.56 

832.0 

0.0000 

-*77.37 

17.24 

0.1868 

7 

0.28040-01 

-0.94650-02 

0. 68230-02 

0.0000 

0.0000 

O.OOOC 

-9 .854 

-0. 88750-01 

8 

0.24690-02 

-0.20750-02 

0.54590-01 

0. 0000 

0. 1117 

0.49761-01 

652.3 

-14.20 

9 

-0.52910-04 

0.70670-03 

-4.575 

-C.2662 

-0. 1 117D-01 

0.1120D-01 

-0. 2464D-01 

0. 2336D-03 

10 

0. 15840-02 

0.81300-02 

4. 374 

0.1331 

0.2234D-02 

-0.9952D-03 

5.040 

O.4630D-O1 

11 

0. 49380-01 

O.OOCO 

0.0000 

-6.655 

-1.117 

0.0000 

157.7 

1.518 

12 

-0.46650-01 

0.22690-02 

-0.10240-01 

-621.3 

-0.5585D-01 

54.75 

1.478 

1.324 

13 

-0. 15770-02 

0.14780-03 

-0. 27290-03 

17.36 

0.0000 

1.000 

0.1035 

C.9342D-03 

14 

-0.19840-04 

0.0000 

-0. 15280-03 - 

-0. 2993D-01 

-3. 088 

0.1119D-02 

0. 3327D-01 

0.0000 

15 

0. 39680-02 

0.0000 

0.0000 

14.98 

-2169. 

0.5598 

-11.09 

-0. 1051 

16 

-0.22040-02 

0.0000 

0.0000 

-8.317 

0.0000 

0.3109 

9185 . 

68.49 
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9 

10 


11 

12 

13 

14 

15 

16 

1 

114.7 

0.0000 


0.65030-01 

-608.1 

-54.62 

9.668 

-0. 2711 

-0. 3949 

2 

-0.25950 

05 -0.12950 

05 

-10.32 

-1.111 

-15.40 

O.OCOO 

0. 0000 

3.896 

3 

4323. 

0.2164D 

05 

14.82 

0.0000 

0.0000 

-27.33 

C. 5109 

-1,116 

4 

81.05 

0.0000 


-0.61260-03 

0. 1592 

-0. 1838 

0.2278 

0.21280-02 

C. 46530-02 

5 

-18.01 

56.07 


0.12250-01 

134.0 

-1871. 

1075. 

22.48 

-0.2791 

6 

-3998. 

37.38 


0.3064 

5.305 

14.71 

-9.111 

0.85150-01 

0.9303 

7 

-3.602 

-1.869 


-0.12250-02 

0.0000 

-1.471 

O.OOCO 

0. 95150-02 

- n. 94 po 

8 

14.41 

3.738 


-0. 19610-01 

-0.4245 

-7.353 

1.822 

0. 1022 

-0. 1116 

9 

-85.73 

1298. 


0.8245 

-0.53060-02 

0.73530-01 

-0.5923 

C.OOOO 

-0.5^820-0? 

10 

0.3602 

-1229. 


-0.8248 

0.42440-02 

-0. 1471 

0. 54670-01 

0.0000 

-0. 29770-02 

11 

-36.02 

8186. 


-14.72 

0.5306 

7.353 

9.111 

C. 85150-01 

0.0000 

12 

- 3.602 

-1.869 


C. 16540-01 

-94.92 

1556. 

1.822 

-18.58 

10. 31 

13 

-0.2161 

0.0000 


0.0000 

•0.1634 

-0.44120-01 

0.91110-01 

-0.9990 

0.9854 

14 

0.0000 

-0.41670- 

•01 

0.27450-04 

0.35820-02 

0.0000 

-54.28 

0. 9989 

0. 83660-03 

15 

0.0000 

-42.05 


0.0000 

2070. 

0.0000 

-123.0 

-n.l oTf 

0.4187 

16 

0.0000 

0.0000 


0.15310-01 

-1169. 

0.0000 

0.0000 

0.0000 

-291. 9 


MATRIX B 

1 

2 

3 

4 

5 

1 

-46.62 

-3.769 

32.91 

-0.32070-02 

0. 18370-02 

2 

-123.5 

-148,0 

-10.82 

-0.63300-02 

-C. 49830-01 

3 

18.83 

1.C65 

0.0000 

0.2811 

0.7788D-02 

4 

-0.1569 

0.1332 

0.0000 

0.0000 

-0.32460-04 

5 

-1.883 

0.0000 

-1406. 

0.90670-02 

0.1456 

6 

-18.82 

2.486 

O.OOCO 

-0.30220-02 

-0.69230-02 

7 

0.8785 

-0.7457 

0.0000 

0. 30220-03 

0.86540-04 

8 

0.2510 

-0.4261 

1.034 

0.18130-02 

0. 17310-03 

9 

0.1883 

-0.17730-02 

-0.77580-01 

0.46850-03 

0.84370-04 

10 

0. 17570-01 

0.6641 

0.62030-01 

-0,24180-04 

0.0000 

11 

0.0000 

0.0000 

O.COOO 

0.0000 

-0.43260-03 

t2 

-0.1882 

0.1065 

0.0000 

-0. 15140-03 

-0.64870-04 

13 

-0.50200-02 

0. 11360-01 

0. 10340-01 

0. 18130- OU 

-0.34620-05 

14 

0.28250-01 

-0. 80110-03 

-60.96 

-0, 13590-04 

C. 28620-03 

15 

-2. 824 

0.0000 

-1454. 

-0.68000-02 

-0.97350-03 

1'&‘ 

O.OCOO 

0.4437 

6.460 

-0.75530-02 

O.OCOO 


MATRIX C 









1 

2 

3 

4 

c 

6 

7 

8 

1 

0.3104D-01 

0.0000 

-0.21840-02 

0.3195 

0.0000 

-0.11940-01 

-0.218 

-C. 3886 0-01 

2 

0. 39500- 05 

0.29050-05 

-0. 15280-03 

-C. 14910-01 

0.0000 

0.0000 

-C. 51240-02 

O.COOO 

3 

0.25460-02 

O.COOO 

-0.2073 

20.26 

523.7 

0.0000 

0, COCO 

-0. *70940-01 

4 

0.0000 

-0.tt153D-03 

10.27 

-2.129 

0.2681 

0.39750-01 

1.^77 

0. 29900-01 

5 

-0. 36680-06 

-0.8298D-07 

-0.12010-04 

-0.21330-03 

-0.62590-04 

O.1196D-0tt 

-0.43910-01 

-G. 20330-0? 

6 

0. 14770-06 

0.0000 

-0. 13010-07 

0.63480-06 

-0.53270-07 

0.23730-07 

-0. 10170-03 

-0.95550-06 

7 

0. 64090-09 

0.42040-08 

O.OOCO 

0.0000 

0.0000 

0.0000 

0,25840-05 

0.24500-0*^ 



9 

10 

11 

12 

13 

14 

15 

16 

1 

0.0000 

O.OCOO 

0.39220-03 

0. 42450-01 

0.0000 

0.1457 

0.0000 

-0.89310-02 

2 

-0.1556 

-0. 41870-01 

0.0000 

-0. 11880-02 

0.0000 

-0.2041D-01 

0.0000 

0.0000 

3 

-51 .78 

53.73 

0.37220-01 

0.0000 

-10.57 

164.7 

-0. 2586 

0.2675 

4 

-5499. 

0.0000 

-0. 39210-02 

0.0000 

-3.530 

5.833 

-0. 13600-01 

0.3870 

c 

0.28820-02 

-0.17930-02 

-0.17650-05 

0. 594CD-04 

0. 35340-03 

0. 16040-02 

-0. 10900-04 

0.33760-02 

6 

0.0000 

0. 17820-05 

0.23370-08 

0.50600-07 

-0.70120-06 

-0. 86890-06 

0. 81210-08 

-0. 17740-07 

7 

0.0000 

-0.26740-04 

-0. 17530-07 

0,0000 

0.0000 

0.0000 

O.OCOO 

O.OCOO 


flATBIX D 

1 2 3 4 5 

1 0,7631 0.0000 -C. 82670-01 0 . 1 45 ID-03 -0 . 4 154D-04 

2 0.0000 1.000 0.1158D-01 -0.6770D-05 0,0000 

3 -5.616 0.0000 2636. 0.434^0-02-0.13140-02 

4 -10.24 0.5673D-01 -1.656 -0. 96690-03 -0.37380-02 

5 0. 10040-03 0.2272D-04 -0.24800-03 0.14500-06 0 . 13860-07 

6 -0.11970-06 -0.67730-07 0.0000 0.28820-09 0.0000 

7 0.0000 0.0000 0.0000 0.0000 O.OCOO 
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Full-State System at Sea- Level Static Condition; Power Lever Angle, 36^ 


MATRIX A 



1 

2 

3 

4 

5 

6 

7 

8 

1 

-2.484 

-0.3262D-01 

-0.2966 

-2556. 

2.969 

604.3 

-1605. 

-13.48 

2 

-0.4666 

-3.034 

196.4 

17.45 

-1.463 

-241.2 

-1339. 

-11. 38 

3 

0.6133D-02 

0.2196D-01 

-90.25 

-25.00 

0.0000 

0.9931 

0.0000 

0.0000 

4 

0. 1124D-02 

0.2013D-02 

0.2026 

-18.33 

0.0000 

-1.076 

1.454 

0. 1207D-01 

5 

-0.3066D- 02 

0.7316D-03 

0.0000 

8. 332 

-249.6 

0.1655 

1.384 

0.4830D-01 

6 

0.3066D-01 

-0.5855D-02 

17.98 

691.6 

0.6991 

-79.12 

24.92 

0.2415 

7 

0.3169D-01 

-0. 1727D-01 

0.0000 

0.8333 

0.0000 

-0.3310D-01 

-17.72 

-0. 1546 

8 

0.1636D-02 

-0.2342D-02 

0.3726D-01 

1.667 

0.2-?96 

0.1986 

7 25.5 

-13.70 

9 

-0.4089D-04 

0.8124D-03 

-4.772 

-0.2500 

-0.3490D-02 

0.1821D-01 

-0.S307D-01 

0. 2415D-03 

10 

0.2887D-02 

0. 1485D-01 

4.556 

0. 5000D-01 

-0.2796D-02 

0.4635D-02 

8.207 

0.6936D-01 

11 

0.5316D-01 

0.0000 

0. 1863 

0.0000 

0.0000 

0.3310 

152.3 

1.352 

12 

-0. 1932D-01 

0.3513D-02 

0.4650D-02 

-529.6 

-0.1049 

57.12 

2.354 

1.168 

13 

-0.6256D-03 

0.2722D-03 

0.5589D-03 

19.63 

-0.1398D-02 

1. 109 

0.1606 

0. 1111D-02 

14 

-0.4599D-05 

0. 9880D-05 

-0. 6286D-03 

0.1876D-01 

-5.359 

0.7447D-03 

0.0000 

-0. 2173D-03 

15 

-0.9202D-02 

-0.3295D-02 

0.2 097 

0.0000 

-2208. 

-0.7447 

-6.230 

-0.5433D-01 

16 

.0000 

0.0000 

0.0000 

10.42 

-0.8738 

0.4138 

8667. 

54.94 



9 

10 

11 

12 

13 

14 

1^ 

16 

1 

-238.9 

196.0 

0. 1131 

-609.7 

-17.89 

-10.98 

0.2190 

-0.6’776 

2 

-0.1839D 05 

-2369. 

-2.651 

0.0000 

17.67 

-21.67 

-0. 4320 

0.4454 

3 

3782. 

0. 1891D 

05 11.11 

-3.931 

-25.30 

O.COOO 

0.3094 

-0.6384 

4 

83.30 

-9.812 

-0.4997D-02 

0.4586 

0.0000 

0.0000 

-0. 7736D-02 

0.0000 

5 

0.0000 

0.0000 

-0. 1 332D-01 

129.7 

-1687. 

1097. 

19.96 

-0.3192 

6 

-3467. 

-184.7 

0. 1466 

15.72 

8.435 

0.0000 

0, 1031 

O.OOCO 

7 

-4.503 

57.72 

0. 3198D-01 

0. 0000 

0.8435 

-1.035 

0. OOCO 

-0. 9789 

6 

0.0000 

9.235 

-0.2665D-02 

0.1310 

0.0000 

6.209 

-0. 2063D-01 

-0. 3405 

9 

-86.00 

1245. 

0.6783 

0.3270D-02 

-0. 421 ID-01 

0.5182D-01 

0. 1546D-02 

0.0000 

10 

0.9006D-01 

-1295. 

-0.7073 

-0.1179D-01 

0.1518 

-0.1242 

0. 1239D-02 

0.34C5D-02 

1 1 

-45.03 

8727. 

-15.03 

O.COOO 

8.434 

-10.35 

0.0000 

0.0000 

12 

0.0000 

-13.85 

0.3330D-02 

-103.8 

1408. 

0.0000 

-16.90 

10. 28 

13 

0.0000 

-0.6926 

-0.5330D-03 

-0.4926 

0.3374D-01 

-0.6209D-01 

-0.9975 

0. 9976 

14 

-0.2027 

0. 1559 

-0.5996D-04 

-0.7371D-02 

-0.5694D-01 

-48.94 

0. 9981 

0. 95'’9D-03 

15 

0,0000 

155.9 

0. 0000 

2040. 

18.99 

93 . 13 

-0.4640 

0.4787 

16 

56,29 

0.0000 

0.0000 

-892.6 

O.OCOO 

-12.93 

0. 1289 

-348.5 


MATRIX B 






1 

2 

3 

4 

5 

1 

-22.73 

-1.507 

-5.475 

-0.2323D-01 

0.9202D-03 

2 

39.58 

-133.8 

10.81 

-0. 1833D-01 

0. 1361D-01 

3 

-3.780 

0.0000 

15,48 

0.2889 

-0. 1300D-02 

4 

0.9452D-01 

0.1242 

-0.1290 

0.0000 

0. 2167D-04 

5 

-0.6299 

-0.1774 

-1030. 

0.2188D-02 

0.1766 

6 

7.560 

2.129 

-5. 161 

0. 0000 

0.2167D-02 

7 

0.3780 

-0.7453 

0.0000 

0.437*70-03 

0.0000 

8 

-0.7560 

-0.2129 

-2.064 

-0. 1751D-02 

-0. 8668D-04 

9 

-0.3151D-01 

0. 1777D-02 

O.OOCO 

0.3720D-03 

-0.2170D-05 

10 

-0. 1260D-01 

0.6182 

-0.5161D-01 

-0. 6127D-04 

0. 4334D-05 

n 

-2.520 

-0.7098 

10.32 

0.0000 

-0. 8667D-03 

12 

-0.1891 

0.2129 

0. 2585 

-0.6562D-03 

-0,433tiD-04 

13 

-0.2772D-01 

0.9938D-02 

0. 1032D-01 

0. 2626D-04 

-0. 3467D-05 

14 

0.8506D-02 

-0.7990D-03 

-51.66 

-0. 9853D-05 

0.2789D-03 

15 

2. 837 

0.7990 

-1103. 

0.9853D-02 

C. 9756D-03 

16 

6.300 

0.0000 

0.0000 

0.0000 

0.2167D-02 


MATRIX C 

1 

2 

3 

4 

5 

6 

7 

8 

1 

0.3670D-01 

-0.2344D-04 

-0.1489D-02 

0.1335 

0.0000 

0.1060D-01 

-9.216 

-0.7689D-01 

2 

0. 1963D-04 

-0.2343D-05 

-0.4471D-03 

0.4000D-01 

0.0000 

0. 1059D-02 

O.OCOO 

O.OCOO 

3 

-0.3271D- 02 

0.0000 

-0.2678 

11.98 

629.2 

-0.9517 

-4.U31 

-0.34720-01 

4 

-0.3274D- 03 

-0.7027D-03 

12.80 

0.0000 

-0.2237 

-0.5302D-01 

0. aa24 

O.OOCO 

5 

0.9827D-07 

-0.2331D-07 

0.893 6D-05 

0.0000 

-0.7834D-04 

-0.2118D-04 

-0.5503D-01 -0.2353D-03 

6 

0. 1924D-06 

-0.4188D-09 

-0.8883D-08 

-0.7947D-06 -0.6667D-C7 

-0. 1263D-06 

-C.13C5D-03 -0. >1100-05 

7 

0. 1365D-08 

0.4886D- C8 

0. 0000 

0.0000 

0,0000 

0.0000 

0.2377D-05 

0. 2073D-07 
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9 

10 

11 

12 

13 

14 

15 

16 

1 

0.0000 

0.0000 

0 . 21340-03 

0.0000 

0.2698 

- 0.4967 

- 0 . 32990-02 

0.0000 

2 

- 0 . 7208 D -01 

0.0000 

0.0000 

0.0000 

0 . 13500-01 

0.0000 

0.0000 

O-OOCO 

3 

64.73 

0.0000 

0.0000 

0.9417 

25.62 

- 46.31 

- 0. 1483 

0.0000 

4 

“ 6131 . 

- 3.698 

0.0000 

- 0. 1049 

- 1.351 

“ 1.657 

- 0 . 49510-01 

- 0 . 34 C 0 D -01 

5 

- 0 . 2160 D -02 

- 0 . 5542 D -02 

0 . 63930-06 

- 0 . 83910-04 

- 0 . 81000-03 

0.0000 

- 0 . 82520-05 

0.4154 0 - 02 

6 

0 . 4294 D -05 

- 0 . 2202 D -05 

- 0 , 25420-08 

0.0000 

- 0 . 24130-05 

- 0 . 19740-05 

0 . 19670-07 

0 . 60880 - 0 ' 

7 

0.00000 

“ 0 . 37430-04 

- 0 . 21610-07 

0.0000 

O.OCOO 

0.0000 

0.0000 

0.0000 


MATRIX D 

1 

2 

3 

4 

5 

1 

0.4032 

0 . 5675 D -02 

0.0000 

0.0000 

- 0 . 69410-05 

2 

0.0000 

1.000 

- 0 . 24770-01 

0 . 14000-04 

- 0 . 13870-05 

3 

0.0000 

- 0.5102 

1556 . 

0.0000 

0.0000 

4 

4.234 

0 . 56700-01 

- 2.478 

- 0 . 14000-02 

0 . 16640-02 

5 

- 0 . 3227 D -03 

- 0 . 56800-04 

0 . 99130-03 

0 . 35020-06 

- 0 . 11100-06 

6 

0 . 3605 D - 06 

0.0000 

- 0 . 14770-05 

- 0 . 25040-08 

0.0000 

7 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 


Seventh -Order System at Sea-Level-Static - Intermediate -Power Condition 






Power Lever Angle, 83® 



MATRIX Ag 








1 

2 

3 

4 

5 

6 

7 

1 

-4.722 

0.1261 

-0.4135 

11.59 

-701. 6 

48C.8 

-527.6 

2 

-0.3438 

-10.65 

191.7 

0.76860-01 

154.9 

-245.8 

-287.0 

3 

0. 59650-01 

1.749 

-35.10 

-0.99370-02 

-27.30 

1.043 

50,84 

4 

-0.7693 

-40.55 

663.9 

-54.74 

423.0 

-125.0 

-747.0 

5 

0.62 080-01 

-0.1598 

1.596 

0.1091 

-1712. 

38.20 

14*’7. 

6 

0.82680-02 

-1.125 

31.36 

-1.489 

5.868 

-52.33 

-10.69 

7 

0.90090-01 

-0.93060-01 

-0.1994 

0.81940-02 

488. 8 

-1.003 

-460.2 

MATRIX Bg 








1 

2 

3 

4 

5 



1 

-75.73 

2.430 

-696.4 

-0.23190-02 

0.33850-02 



2 

5. 136 

-502.9 

27.63 

0.21960-01 

0.41230-02 



3 

-1.538 

75.31 

-22.72 

0.2501 

-0. 22480-02 



4 

-36. 16 

-1697. 

-362.9 

0.1751 

0.20970-01 



5 

12.35 

2.814 

-3731. 

-0. 10490-01 

0.61100-01 



6 

-0.20750-01 

-43. 14 

-13.50 

0. 37530-02 

-0.67670-03 



7 

1.441 

-5.879 

501.4 

0.96860-03 

-0.50710-03 



MATRIX Cg 








1 

2 

3 

4 

5 

6 

7 

1 

0. 18590-01 

-0.10650-03 

0.66090-03 

0. 35320-03 

-0. 1099 

-0.10590-01 

0.1167 

2 

-0.1122D-06 

0.91370-05 

-0.20060-03 

0.21010-04 

0.12010-01 

-0. 57460-03 

-0. 12740^ 

3 

-0. 18960-01 

-0.42890-01 

0.4850 

-0.3405 

454.9 

-1.552 

100.1 

4 

-0. 29190-01 

-1.075 

17.19 

0. 15880-01 

26.33 

-0.8305 

-36.16 

5 

0.24160-05 

-0.16360-04 

0. 15480-03 

-0.25620-05 

-0.1538 

0.32280-03 

0. 1414 

6 

0. 13370-06 

-0.90950-11 

-0. 12130-07 

-0. 66530-11 

-0. 41860-07 

0.13180-09 

-0. 15090- 

7 

0.56190-09 

0.37500-08 

-0.63890-07 

-0. 73400-13 

0.65930-08 

0. 33060-10 

0. 18420- 

MATRIX Dg 








1 

2 

3 

4 

5 



1 

0.6905 

-0.67110-01 

-0.38580-01 

-0.97370-06 

0.41830-04 



2 

-0.10200-03 

1,0000 

0.13320-01 

0.96700-07 

-0.28740-05 



5 

2.858 

-3.175 

5766. 

0. 14810-C2 

-0. 20610-01 



4 

-0.2215 

-47. U2 

11.46 

0. 49240-02 

0.75110-03 



5 

0.12250-02 

0.26340-03 

-0.1595 

-C. 11010-05 

0.37520-07 



6 

0.60050-10 

-0.84900-09 

-0.43540-07 

0.285SD-11 

0.43530-12 



7 

-0. 12010-09 

-0.14460-07 

0.44840-08 

-0.14300-12 

-0.89340-13 




30 


MATEIX Hg 



1 

2 

3 4 

5 

6 

7 

1 

1.000 

0.0000 

0.0000 0.0000 

0.0000 

0.0000 

0, 0000 

2 

0.0000 

1.000 

0.0000 0.0000 

0.0000 

0.0000 

0.0000 

3 

0.0000 

0.0000 

1.000 0.0000 

0.0000 

0.0000 

0. 0000 

4 

0.0000 

0.0000 

0.0000 1.000 

0-0000 

O.OCOO 

0.0000 

5 

0.0000 

0.0000 

0.0000 0.0000 

1.000 

0.0000 

O.OOCO 

6 

0.0000 

0.0000 

0.0000 0.0000 

0-0000 

1,000 

0.0000 

7 

0.0000 

0.0000 

0.0000 0.0000 

0.0000 

0.0000 

1.C00 

8 

0.9825D-04 

0.2172D-03 

-0. 1596D-02 0. 1589D-03 

0.9998 

-0-3147D-02 

7.006 

9 

-0.2913D-01 

-1.082 

17.30 0.1446D-01 

29.21 

-0,7879 

-38.97 

10 

0. 1733D-03 

0.4969D-02 

0.9634 0.1744D-05 

-0. 9705D-01 

-0.5375D-03 
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table I. - EIGENVALUES FOR FULL STATE 
LINEAR MODELS AT SEA LEVEL STATIC 


AND VARIOUS POWER CONDITIONS 


Power 
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angle, 

deg 

Eigenvalue, 

rad/sec 
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angle, 

deg 
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-51.24±j312.73 
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-292.9 
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-1344.9 


TABLE n. - EIGENVALUES FOR FULL-STATE AND 
REDUCED-ORDER LINEAR MODELS AT SEA LEVEL 
STATIC INTERMEDIATE POWER CONDITION 


Power lever angle, 83°. 


Linear 

model 

Eigenvalue, 

rad/sec 

Linear 

model 

Eigenvalue, 

rad/sec 

Full state 
(16th order) 

-4.86 

-6.49 

-15.22±jl0.75 
-20.74 
-37.02±]8.79 
-89. 3±j40. 36 
-I25±j55. 31 
-49. 4±jl24. 8 
-38. 1±J333. 9 
-1542.6 

7th order 

-49.3 
-73. 1 
-2141 

5th order 

-4.25 

-6.85 

-38.84±jl6.95 

-154.4 

3rd order 

-4. 24 
-6.61 
-39.8 

7th order 

-4.80 
-7. 52 
-22.4 
-31. 1 


High- pLow-pressure 



Rgure 1. - Schematic representation of noo-PW-100t3) augmented turbofan engine. 
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Figure 2 . - Computational flow diagram of real-time FlOO-PW-100 engine simulation. 



Figure 3. - Generation of data for linear models. 
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Figure 4 - Sampling and recording processes for state variable disturbance data (SEL810B computer). 



Figure 5. - Data verification, processing, and linear mode! generation. 
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Figure 6. - Comparison of transient response of nonlinear hyorld simulation and full-state model for step change In exhaust- 
nozzle area. Sea-level -static - Intermediate-power condition. 



Change in exhaust nozzle pressure, P7. N/cm^ Change in burner pressure, P4, N/cm^ Change in fan speed, N|^, rpm 







(e) Exhaust-nozzle pressure. (f) Thrust, 

Figure 8. - Comparison of linear models for step change in exhaust nozzle area. Sea- level -static 
intermediate-power condition. 
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Change in exhaust nozzle pressure, Py, N/cm^ Change in burner pressure, P 4 , N/cm'^ Change in fan speed, N|_, rpm 


80 ,- 



(e) Exhaust nozzle pressure. (f) Thrust. 

Figure 9. ~ Comparison of linear models for step change in main burner fuel flow. Sea-level-static - 
intermediate-power condition. 
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Figure 10. - Eigenvalues for fulhstateand fifth-order linear models, 
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